On the formation of non-topological string networks 



Ana Achucarro 
Department of Theoretical Physics, UPV-EHU, Bilbao, Spain 
Institute for Theoretical Physics, University of Groningen, The Netherlands 

Julian Borrill 

Center for Particle Astrophysics, University of California, Berkeley, CA 94720 
National Energy Research Scientific Computing Center, Lawrence Berkeley National Laboratory, University of California, 

Berkeley, CA 94720 

Andrew R. Liddle 

Astronomy Centre, University of Sussex, Falmer, Brighton BN1 9QJ, United Kingdom 

(February 1, 2008) 

We review recent simulations of the formation of a particular class of non-topological defects known 
as semilocal strings during a phase transition. Semilocal strings have properties that are intermediate 
between topological cosmic strings and stable electroweak strings, and therefore the observation that 
they could form in substantial quantities during a cosmological phase transition has implications 
for structure formation, baryogenesis, etc. At the same time, and from a purely theoretical point 
of view, they provide a very good testing ground for investigating the role of gauge fields in defect 
formation. 
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I. INTRODUCTION 

Our current understanding of particle physics and cos- 
mology implies that the early Universe probably went 
through a series of symmetry breaking phase transitions 
as it cooled down and expanded to become what we 
know today. In these transitions topological (and possi- 
bly non-topological) defects are expected to have formed 
Although the cosmological evidence for the exis- 
tence of such defects remains unclear there is plenty 
of experimental evidence from condensed matter systems 
that networks of defects do form in symmetry breaking 
phase transitions according to the Kibble-Zurek mecha- 
nism 0,^,^], the most recent confirmation coming from 
the experiments in vortex formation in superfluid Helium 
that have been performed here in Helsinki at the ULTI 
and in Grenoble and Lancaster |6|. 

From the theoretical point of view, the last few years 
have seen great progress in our understanding of the cre- 
ation of topological defects in a symmetry breaking phase 
transition, largely through the interaction between cos- 
mologists and condensed matter physicists in a beautiful 
example of interdisciplinary collaboration ||. However, 
it is fair to say that so far most of this progress has been 
in cases involving only scalar fields and global symme- 
tries, whereas most theories of elementary particles are 
based on gauge, rather than global, symmetries, and it 
is these models that are relevant to cosmology. 

Defect formation in these systems remains an open 
problem. For a start, the usual order parameter (the 
expectation value of the Higgs) is not gauge-invariant, 



and it is not obvious in general how to choose a gauge- 
invariant order parameter. The applicability of the 
geodesic rule is also under scrutiny . At a more basic 
level, even the nature of the electroweak phase transition 
is under discussion. Under the circumstances, it is dif- 
ficult to make reliable estimates of the initial density of 
gauged defects such as electroweak || or semilocal [^0 11 
strings, which is of course a necessary first step towards 
understanding their cosmological implications for struc- 
ture formation, baryogenesis or generation of primordial 
magnetic fields. 

Gauge fields are often neglected in analytic calcula- 
tions, largely because the usual techniques do not lend 
themselves easily to their inclusion. But if we are to 
believe numerical simulations we have to conclude that 
this may not be such a good idea; whenever present, 
gauge fields tend to play a crucial role in the dynam- 
ics. A rather drastic example is the existence of non- 
homogeneous ground states (even in the absence of ex- 
ternal magnetic fields) in large samples of certain two- 
dimensional unconventional superconductors whose or- 
der parameter has a magnetic moment. An incommen- 
surate phase appears which can be described as a paral- 
lel array of (non-topological) domain walls; the magnetic 
field lives in the core of these walls, where there is sym- 
metry restoration, with the magnetic flux alternating in 
sign from one wall to the next so that the net magnetic 
flux through the sample remains zero. These inhomo- 
geneous phases have not been observed experimentally, 
as there is no known system in the right range of values 
of the Landau-Ginzburg parameters, but the theoretical 



1 



prediction is well established 

Another example of the effect of gauge fields on the 
dynamics is the observed breakdown of the geodesic rule 
(the assumption that the field configuration will take 
the shortest interpolating path between different vacuum 
states) during first order phase transitions in systems 
with gauge invariance [____]. A recent numerical study 
by Copeland and Saffin, both in the case of the abelian 
Higgs model and the 9 W = limit of the Wcinbcrg-Salam 
model, has found there may be enough energy in the col- 
lisions of separated nucleated bubbles to drive the Higgs 
value away from its vacuum expectation value and make 
it oscillate around zero, so extra regions with symmetry 
restoration are produced. The validity of the geodesic 
rule seems to depend not only on the initial bubble sep- 
aration but also on the gauge coupling, indicating gauge 
fields are playing a crucial role. 

The upshot of both studies would seem to be that the 
backreaction of the gauge fields on the scalars can re- 
sult in an enhancement of symmetry restoration leading 
to larger defect densities than one would have inferred 
from the scalar fields alone. Here we present recent re- 
sults (l_]|l_| that corroborate this idea in the context of 
electroweak strings. 

In particular, we have studied the formation of a type 
of non-topological vortices known as semilocal strings; 
these are stable Abrikosov-Niclscn-Olesen vortices in the 
9 W — > 7r/2 limit of the (bosonic sector of the) Weinberg- 
Salam model pj. Note that in these models the vacuum 
manifold -the manifold of ground states- is simply con- 
nected thus if one ignored the effect of gauge fields one 
would not expect any vortices to be formed at all. 

Before we continue, we should point out that analogous 
systems have been studied in condensed matter. In |L5"| |, 
the system was an unconventional superconductor where 
the global SU(2) group was the spin rotation group. In 
[M the hypothetical case of an "electrically charged" A- 
phase of 3He (a superconductor with the properties of 
3He-A) was considered. In this case the global group 
was SO (3), the group of orbital rotations. Both papers 
discussed textures (continuous vortices) in such super- 
conductors, which correspond to the "skyrmionic" con- 
figurations of |TT|j2^ 1 (see Fig. 1). 

We have carried out a numerical study of semilocal 
string formation through the Kibble-Zurek mechanism 
based on a novel numerical technique |l3|Jl4| which can be 
applied to other gauged non-topological defects as well. 
The strategy we propose is to follow the evolution of the 
field strength, which is used as an indicator for the pres- 
ence of defects; the initial conditions are obtained by an 
extension of the Vachaspati-Vilenkin algorithm jlTj ap- 
propriate to non-topological defects, plus a short period 
of dynamical evolution including a dissipation term (nu- 
merical viscosity) to aid the relaxation of configurations 
in the 'basin of attraction' of the defects under study. 

We observe that, even if we assume no symmetry 



restoration in the initial stages of the phase transition, 
the interaction between the gauge fields and the scalar 
fields is such that symmetry restoration will eventually 
occur simply because it is energetically favourable for 
the magnetic field to be concentrated in regions where 
the Higgs field has a value close to that of the symmetric 
phase. Thus, vortices are formed. 

Since we are proposing a new technique, we test it 
in several ways. Restriction to the abelian Higgs model 
gives excellent agreement with analytic and numerical 
estimates for cosmic strings, a simpler type of topolog- 
ical defect analogous to a superconduting flux-tube, as 
given in (l7), and we also test the robustness of our re- 
sults under varying initial conditions and numerical vis- 
cosity. Finally, our results seem to be in good agreement 
with previous estimates for semilocal string formation in 

mi. 

It is worth mentioning that our analysis suggests 
a similar behaviour for electroweak strings away from 
9 W = it/2, as long as they are stable. Nagasawa and 
Yokoyama p^j estimated the initial density of (unstable) 
electroweak strings by looking at the thermal distribution 
of scalar fields and concluded that the density would be 
negligible. While their calculation at sin 2 9 W ~ 0.23 is 
not necessarily in contradiction with ours at 6 W = 7r/2, 
their method neglects the role of the gauge fields, which 
we would argue to be a key ingredient. Thus, full numer- 
ical simulation appears to us to be the only reliable way 
to investigate electroweak string formation rates. 

In what follows we briefly summarize our tech- 
nique and results. For more information we refer 
the reader to our papers |T^,|l4],|l| and web pages. 
Colour images and movies of the three-dimensional 
si mulations can currently be found on the WWW 



at tittp: / / cfpa.berkeley.edu / ^borrill / defects / semilocal.html 
Also, a 100 frame mpeg movie (0.5Mb) of the two- 
dimensional simulation can be viewed at 



http://star-www.cpes.susx.ac.uk/people/arl_recent.htm' 



II. THE MODEL 

We work in flat space-time throughout. The semilocal 
model is described by the following lagrangian: 

C = (V M - ieA^) 4 (V + i&A") fa 
+ (V„ - *eA M ) 0_ (V+teA")0_ 
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(i) 



where fa and fa are two equally-charged complex scalar 
fields, is a U(l) gauge field and F^ v = \7^A„ — V^A^ 
the associated gauge field strength. Notice that this is 
just the bosonic sector of the Weinberg-Salam model in 
the 9 W = ir/2 limit in which the SU(2) coupling constant 
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is set to zero and the W-bosons decouple. In particu- 
lar its vacuum manifold is the three-sphere S 3 , which 
has no non-contractible loops, and yet stable Abrikosov- 
Nielsen-Olesen vortices fl22j23f can form @Q. On 
the other hand, setting 02 = obtains the abelian 
Higgs model; thus, comparison with topological strings 
is straightforward, and we will use it repeatedly as a test 
case, both to check our simulation techniques and to min- 
imize systematic errors when quoting formation rates. 

As in the abelian Higgs model, the gauge coupling 
and the vacuum expectation value of the Higgs can be 
rescaled to one by choosing appropriate units (the in- 
verse vector mass, for length, and the symmetry breaking 
scale, for energy) . In that case the Lagrangian becomes 



+ (V M - iA, 



- U\tl\ 2 + \fo\ 2 - 



(2) 



and the only remaining parameter in the theory is j3 = 
mg/rriy, the ratio between the scalar and vector masses 
(squared), whose value determines the stability of an in- 
finitely long, straight, string with a Nielsen-Olesen pro- 
file: it is stable for j3 < 1, neutrally stable for j3 = 1 and 
unstable for f3 > 1 |lO|JHlJ]. For (3 = 1 there is a fam- 
ily of solutions with the same energy and different core 
widths, of which only the semilocal string has complete 
symmetry restoration in the center [Tl|p^ |. 

Whenever (3 < 1 we expect a certain amount of these 
non-topological strings to form during a phase transi- 
tion; the lower the value of (3, the higher the formation 
rate. The enhancement of symmetry restoration is eas- 
ily understood from the minimal coupling between the 
gauge field and the scalars. In regions where the mag- 
netic field is strong (because of fluctuations), the term 
A 2 (j) 2 in the lagrangian will tend to drive the Higgs' ex- 
pectation value towards zero, in competition with the 
potential term. This is basically the same physics as 
in the Meissner effect, the expulsion of magnetic fields 
from a superconductor. In finite samples subjected to 
an external magnetic field the system reacts by forming 
surface currents that screen the magnetic field, but if 
the magnetic field is strong enough it will penetrate the 
sample through regions where the symmetry is restored 
(e.g. Abrikosov vortices, in type II superconductors [B2fl). 
Similarly, in a cosmological context, the best way to cope 
with a strong fluctuation of the magnetic field is to reduce 
the value of $ towards zero. 



~ VVa +^W + ^+^+^J- 1) 1>a (3) 

+A 2 ip a + {-l) b (2A ■ V + V • A) ip b = , 

(where b is the complement of a — I <-> 2, 3 <-> 4 and 
dots are time derivatives) for the scalar fields and 

At - V 2 Ai + V 4 V • A + 2 ( tpi V« ip 2 + V> 3 Vi ^4 ) (4) 



+2A, +i> 2 2 + V> 3 2 + i/>l) = , 

for the gauge fields (i = 1,2,3), together with Gauss' 
law, which here is a constraint derived from the gauge 
choice, and is used to test the stability of the code, 



2 Ui V + Vs V ip* + V* A* = . 



(5) 



The arrows indicate asymmetric derivatives. 

This system is discretized using a standard nuerical 
technique known as the staggered leapfrog method. How- 
ever, to reduce its relaxation time we also add an ad hoc 
dissipation term to each equation (j]i^ a and rjAi respec- 
tively). In an expanding Universe the expansion rate 
would play such a role, though 77 would typically not be 
constant. We tested a range of strengths of dissipation, 
and checked that it did not significantly affect the num- 
ber densities obtained. The simulations we display later 
used 77 = 0.5 and periodic boundary conditions. 

We estimate the number density of defects by an ex- 
tension of the Vachaspati-Vilenkin algorithm ]l7|] : we 
first generate a random initial configuration for the scalar 
fields drawn from the vacuum manifold, which is not dis- 
cretized, then find the gauge field configuration that min- 
imizes the energy associated with (covariant) gradients^. 
If space is a grid of dimension N x N x N, the cor- 
relation length is chosen to be some number p of grid 
points (p = 16 in our simulations). To obtain a reason- 
ably smooth configuration for the scalar fields, we throw 



down random vacuum values on a 



N 



N 



subgrid; 
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the scalar field is then interpolated onto the full grid by 
bisection. Strings are always identified with the location 
of magnetic flux tubes. For cosmic strings this accurately 
reproduces the standard results Jf7| . 

For semilocal strings, on the other hand, the initial 
configurations thus generated have a complicated flux 
structure with extrema of different values (top panel of 
Fig. 1), and it is far from clear which of these, if any, 
might evolve to form Nielsen-Olesen vortices; we resolve 



III. NUMERICAL SIMULATIONS 

We work in temporal gauge A§ = 0. Splitting up the 
scalar fields into four real scalars via <f>\ = %p\ + it^i-, 
4>2 = ip3 + the equations of motion are 



*In fact, the energy-minimization condition is redundant, 
since the early stages of dynamical evolution carry out this 
role anyway. Having checked that our results were indepen- 
dent of the initial conditions for the gauge fields, we used 
Ai(x) = ij>iViip2 — faViipi + V>3Vi^/>4 — ip4Viil>3 as the initial 
gauge field configuration in our simulations. 
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this ambiguity by numerically evolving the configurations 
forward in time. As anticipated, in the unstable regime 
ft > 1 the flux quickly dissipates leaving no strings. By 
contrast, in the stable regime ft < 1 stringlike features 
emerge when configurations in the "basin of attraction" 
of the semilocal string relax unambiguously into vortices 
(bottom panel of Fig. 1). 

The correlation length in the simulations is only con- 
strained to be larger than the size of the vortex cores, 
to avoid overlaps. This results in a minimal value of the 
parameter ft of around 0.05 ^. For each of seven different 
values of ft, we take several initial configurations on a 64 3 
(or 128 2 ) grid smoothed over every 16 grid-points. As ex- 
pected, for ft < 1 we find a formation rate which depends 
on ft, tending to zero as ft tends to 1. The formation rate 
is lower than in the cosmic string case, but significant: 
at the smallest ft we were able to simulate, ft — 0.05, the 
rate is a fraction of order 1/3 of that of cosmic strings. 
Our statistical results are derived from a large suite of 
simulations (700 in all) carried out on a 64 3 grid using 
a SUN Ultra II workstation. We have also performed a 
large 256 3 simulation on the Cray T3E at NERSC from 
which we have extracted the image in Fig. 4. 



IV. RESULTS 

Before discussing the results in any more detail, a 
few comments are in order about the general procedure 
we are following. We checked our strategy on a two- 
dimensional toy model obtained by ignoring one of the 
spatial coordinates Jl3[ |. 

We first tested the performance of our codes for the 
case of cosmic strings by ignoring one of the (complex) 
scalar fields, setting ip3 = tp^ = 0. This makes the defects 
topological, and the flux tubes formed now map out the 
locations of winding in the scalar field configurations. We 
obtained a cosmic string number density 



= 0.32 ±0.02 



(6) 



per correlation area, in perfect agreement with the ana- 
lytic value of 1 /3 obtained with the standard Vachaspati- 
Vilenkin algorithm on a 2-dimensional square lattice 
without discretizing the vacuum manifold. The good 
agreement arises because the simulations are themselves 
carried out on a square lattice, and confirms that our 
identification of cosmic strings via the flux tube struc- 
ture works extremely well [[[ 




FIG. 2. The number of semilocal strings formed per initial 
two-dimensional correlation volume. Each point is an average 
over ten simulations. Squares indicate that the vacuum initial 
conditions described in the text were used, while open circles 
indicate that non-vacuum (thermal) initial conditions were 
used. 



We also tested the robustness of our conclusions under 
varying initial conditions. We have already mentioned 
that the initial configuration for the gauge fields was 
seen to play no appreciable role. On the other hand, 
in a second-order phase transition one expects the scalar 
field to be out of the vacuum in some sort of thermal 
distribution. This has been studied by Ye and Branden- 
bcrgcr [^7) , and we closely followed their strategy, draw- 
ing the scalar field magnitude randomly from a gaussian 
distribution instead of the uniform distribution they use. 
Fig. 2 shows the results of two-dimensional simulations 
with two sets of initial conditions, one with no symmetry 
restoration and other with a more thermal distribution 
of scalar field values. We examined different choices for 
the width of the gaussian, and concluded that any rea- 
sonable choice makes no difference to the results; in those 
displayed in Fig. 2, the gaussian has a dispersion equal 
to the vacuum expectation value of the scalar field. 

As can be seen, the influence of initial conditions is 
minimal. The same general shape for the /3-dependence is 
found, and the number of strings identified per simulation 
is extremely close, with the error bars overlapping. 

Our results are also in good agreement with those of 
, in which simple configurations some way away from 
the semilocal string were permitted to relax into semilo- 
cal strings. Finally, we note also that extrapolation of 
Fig. 2 to ft — is consistent with the estimated value of 



T If ft is lowered further, the scalar string cores become too 
wide to fit into a correlation volume, in contradiction with the 
vacuum values assumed in a Vachaspati-Vilenkin algorithm. 

*The error quoted is the standard error on the mean, but 



note that the uncertainty from discretization is somewhat 
larger: on a triangular lattice the Vachaspati-Vilenkin algo- 
rithm would give a formation rate of 1/4 per triangle rather 
than 1/3 per square J26j , though one would also have to ac- 
count for the difference in area. 
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FIG. 1. The flux tube structure in a two-dimensional semilocal string simulation with f3 = 0.05. The upper panel (t — 0) 
shows the initial condition after the process described in the text. The lower panel shows the configuration resolved into five 
flux tubes by a short period of dynamical evolution (t = 100). These flux tubes are semilocal vortices. Note the different 
numbers of upward and downward pointing flux tubes, despite the zero net flux boundary condition. The missing flux resides 
in the smaller 'nodules', made long-lived by the numerical viscosity; the expansion of the universe could have a similar effect 
and preserve these 'skyrmionic' configurations 
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1/8 w 0.12 found in Ref. g^. 

This is a two-dimensional simulation, but it is also ex- 
pected to be representative of any cross section of the 
three-dimensional one. 

We can now proceed to describe the results of our sim- 
ulation in some detail. 

In order to bring down the statistical errors on our 
result we carried out a large number of runs; the counting 
of flux tubes was automated by adopting the criterion 
that any isolated extrcmum in the magnetic flux which 
was at least half that of a perfect Nielsen-Olesen vortex 
was identified as a string. We tested this criterion on 
selected simulations which we studied in detail ourselves, 
and concluded it gave an accurate counting. 

As in the two-dimensional simulations (Fig. 1), the 
starting configurations obtained by our described proce- 
dure initially yield a complicated mess of flux. However, 
after a few timesteps the flux resolves itself into loops and 
open segments of string. We observed a clear interaction 
between nearby segments which join to form longer seg- 
ments, as conjectured in |24|,[]iJ. Fig. 3 shows an open 
segment of string whose ends meet to form a loop. 

Fig. 4 shows two time slices from a single large (3 = 0.05 
simulation using the Cray T3E at NERSC; these are 
zooms showing only part of the 256 3 simulation box. We 
see a collection of short string segments and loops; visu- 
ally this is very different from a cosmic string simulation 
where strings cannot have ends. As time progresses, the 
short segments either disappear or link up to form longer 
ones. 

We can immediately conclude from these images that 
the formation rate of semilocal strings is not extremely 
close to zero; the fact that flux tubes are observed in 
our simulations implies that the formation rate cannot 
be much smaller than one per correlation volume. 

In order to quantify the formation rate, we compute 
the total length of string in the simulations, always com- 
paring the semilocal string density to that of a cosmic 
string simulation with the same properties (including dis- 
sipation) to minimize systematic errors. We determine 
the length by setting a magnetic flux threshold and com- 
puting the fractional volume of the box which exceeds 
it. In Fig. 5, we plot the length of semilocal string rel- 
ative to the length found in cosmic string simulations, 
as a function of time and with (3 = 0.05. We see that 
after a transient during which the initial tangle of flux 
sorts itself out, the system settles down to a reasonable 
equilibrium. The upward trend appears to be caused by 
the periodicity of the simulation box, which freezes-in 
any string crossing the box, favouring cosmic string an- 
nihilation because of their higher density. We take the 
relative densities of semilocal and cosmic strings to be 
that at time 50 in these simulations. There is a modest 
dependence on the choice of flux threshold, and we set it 
at one-half the flux density of a Nielsen-Olesen vortex. 

Fig. 6 shows the resulting ratio of semilocal and cos- 




FIG. 3. Two snapshots, at t = 70 and t = 80, of a 64 3 
simulation with = 0.05 where the ends of an open segment 
of string join up to form a closed loop. The loops in our 
simulation seem to behave like those of topological cosmic 
string, contracting and disappearing. 

mic string lengths, as a function of the stability param- 
eter (3. These results are derived from 50 simulations 
(both semilocal and cosmic) at each /3 value, carried out 
in boxes of dimension 64 3 . The error bars include the 
statistical spread between simulations, and an estimated 
25% systematic error from the length counting algorithm 
(see the spread in Fig. 5) and the viscosity. Those lat- 
ter uncertainties are the dominant ones. Recalling that 
the formation rate of cosmic strings is estimated to be 
of order one per correlation volume (0.88 in Ref. 0), 
these results are in excellent agreement with the two- 
dimensional results we found in Fig. 2. They show a 
significant formation rate for low /?, decreasing dramat- 
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FIG. 4. Part of the large simulation, shown at time 60 and 
time 70. Note several joinings of string segments, e.g. two 
separate joinings on the long central string, and the disap- 
pearance of some loops. The different apparent thickness of 
strings is entirely an effect of perspective. 



ically as (i — > 1, beyond which there is no evidence of 
semilocal string formation. The small amount of string 
seen in some large-/? simulations is an artifact of the vis- 
cosity; the flux all dissipates if the viscosity is turned off, 
while it persists if /3 < 1. 



V. CONCLUSIONS 

The formation of defects during a symmetry-breaking 
phase transition in systems with gauge invariance is a 
very important issue in cosmology, but it is not as well 
understood as it is in systems with only global symme- 
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FIG. 5. This shows the ratio of total string lengths in a 
semilocal and cosmic string simulation, with f3 = 0.05. The 
different lines show different magnetic flux thresholds, from 
bottom to top they are 0.6, 0.55, 0.5, 0.45 and 0.4 times the 
peak flux of a Nielsen-Olesen vortex. 
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FIG. 6. The ratio of lengths of semilocal and cosmic 
strings. 



tries. The lack of experimental data from condensed mat- 
ter systems forces us to rely on numerical simulations to 
test our ideas. We have presented here a study of non- 
topological defect formation where the role of gauge fields 
is seen to be crucial and leads to a very different scenario 
from what one would have inferred from the scalar fields 
alone. 

Our simulations have shown that dynamically sta- 
ble non-topological defects, such as semilocal strings, 
could form in substantial quantities during a cosmological 
phase transition. This seems to happen even if the con- 
figuration immediately after the phase transition has no 
appreciable symmetry restoration (recall our initial con- 
ditions place the scalar field in the vacuum everywhere), 
and is due to the back-reaction of the gauge fields on the 
scalars, which favours a low vacuum expectation value 
for the latter. In the case of non-topological defects, the 
intensity of the effect, and therefore the number density 
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of defects created, increases with the classical stability of 
the configurations. For semilocal strings with (3 = 0.05 
we found a formation rate as high as one third of that 
of topological strings, while as (3 was increased towards 
the stability/instability transition at (3 — 1, the density 
dropped to zero. 

Moreover we have observed short segments of string 
growing and joining to form longer strings and loops, as 
conjectured in |l^,|2^], while others collapsed longitudi- 
nally, as invoked in some baryogenesis scenarios ]2S| ] . The 
details of these dynamics seem to depend on the length of 
the segments and the distance between neighbouring ones 
and deserve further study, as do the scaling behaviour of 
the network and the implications for structure formation, 
baryogenesis, etc. 

We suspect that a similar behaviour is to be expected 
of (stable) electroweak strings. Nagasawa and Yokoyama 
Ipof studied electroweak string formation and concluded 
that the initial density would be totally negligible for re- 
alistic values of the weak mixing angle. However, this is 
not necessarily in contradiction with our results for two 
reasons. First of all, in our simulation semilocal strings 
arise during the evolution because of the back-reaction 
of the gauge fields on the scalars, which enables initially 
short pieces of string to grow and join up to form longer 
ones, an effect not included in their analysis. Secondly, 
their calculation takes place at sin 2 9 W ~ 0.23, for which 
the electroweak string is dynamically unstable p9j; in 
this case we would find that all the flux dissipates soon 
after the phase transition. In any case, numerical simu- 
lation appears to be the only reliable way to address this 
problem. 

In the meantime, the semilocal string model remains 
an excellent testing ground for defect formation scenarios 
in the presence of gauge symmetries. 
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